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Abstract 

We establish exactly solvable models for the motion of neutral parti- 
cles, electrically charged point and spin particles (U(l) symmetry), isospin 
particles (SU(2) symmetry), and particles with color charges (SU(3) sym- 
metry) in a gravitational wave background. Special attention is devoted 
to parametric effects induced by the gravitational field. In particular, 
we discuss parametric instabilities of the particle motion and parametric 
oscillations of the vectors of spin, isospin, and color charge. 



1 Introduction 

Periodic external fields are known to induce parametric phenomena in physical 
systems. This includes such effects as parametric oscillations (the oscillation 
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frequency becomes a periodic function of time) and parametric instabilities (ex- 
ponential growth of certain dynamical quantities) El ■ Classical examples 
are the parametric resonance in vibrations of mechanical and electrical sys- 
tems m 121 E] and plasma instabilities in external electromagnetic fields ^ E] . 
From a mathematical point of view, such phenomena are described by differen- 
tial equations with periodic coefhcients which are subject to Floquet's theory 
m El El 13 m El ■ As we shall demonstrate here, equations of this type naturally 
appear if one considers the motion of different kinds of particles in a gravita- 
tional wave (GW) background. The gravitational wave may play the role of 
an external periodic pumping field. In fact, most of the attempts for a direct 
detection of gravitational waves rely on this concept. The mathematical similar- 
ity to dynamical equations which are known to describe parametric oscillations 
and resonances then naturally suggests the possibility of gravitationally induced 
parametric effects. 

The idea of parametric phenomena in GW fields is not new. For a linearized 
GW field it was discussed in JD). The first exactly solvable model for the 
evolution of a kinetic system in a nonlinear GW field, demonstrating explicitly 
the possibility of parametric excitation of a relativistic plasma by a periodic 
GW, has been established in ^Jj. At the same time, in the 90th, the problem 
of parametric resonance during the reheating phase of inflationary models has 
become an intensely elaborated topic in a cosmological context (see, e.g., [T^IT^ 

15, I6|). More recently, investigations, concerning parametric phenomena in 
a GW field, have attracted attention again (see, e.g., jElElEI)- 

The purpose of the present paper is to clarify characteristic features of grav- 
itationally induced parametric effects for simple dynamical configurations. As 
a first example we consider the motion of an electrically charged point particle 
which is simultaneously exposed to a constant magnetic field and a gravitational 
wave with front plane orthogonal to the magnetic field. In the second exam- 
ple we include an additional spin degree of freedom which is described by the 
Bargmann-Michel-Telegdi (BMT) equations [201 • In the third case the electri- 
cally charged spin-particle is replaced by a particle with isospin and the magnetic 
field is replaced by a corresponding Yang-Mills field. The isospin dynamics is 
governed by Wong's equation j^J for the three-dimensional isospin vector. Fi- 
nally, we consider the motion of particles with color charge, described by Wong's 
equation for the eight color degrees of freedom. Using suitably specified Yang- 
Mills fields, we establish a general scheme which allows a unified treatment of the 
particle dynamics for all four cases. We show that on a periodic GW background 
this dynamics is characterized by Hill and Mathieu equations. Well-known sta- 
bility properties of the latter allow us to classify the particle motion accordingly. 
This implies parametric oscillations and/or parametric instabilities as generic 
phenomena. The precession dynamics of the vectors of spin, isospin and color 
charge is coupled to the particle motion and parametrically driven as well. 

The paper is organized as follows. In Section El we establish the basic dy- 
namic equations for the Abelian and non-Abelian subcases to be discussed in 
the following. In SectionElthe Yang-Mills fields for the latter cases are specified. 
The (exact) gravitational background is characterized in Section 01 Section El 
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is devoted to a compact, general solution for the particle dynamics. A "sand- 
wich" GW is considered as a special case. The spin precession for an electrically 
charged particle is the subject of Section Sections [T] and |S1 discuss the dy- 
namics in the spaces of isospin and color charge, respectively. In Section |3 we 
summarize our main results. We use units in which h = c = 1. 



2 Particle dynamics: Basic equations 

Let us consider the evolution of relativistic point particles with either an electric 
charge and a spin-vector, or an isospin, or a color charge. The concepts of 
classical particles with isospin (for the SU(2) symmetry) or color charge (for 
the SU(3) symmetry) are generalizations of the electrically charged particles to 
the non-Abelian case (see, e.g., ED). The dynamical equations for the particle 
momentum p^, for the spin- vector S"* and for the charge Q^'^\ where (A) is a 
group index, are 

^^g^ and ^ = gi-K (1) 

respectively. Here, D denotes the covariant differential, r is a parameter along 
the particle worldline and J-"' is the force four vector which is orthogonal to the 
particle momentum = mdx''/dT, i.e., piT"^ = 0. The quantity describes 
the spin rotation within the BMT theory. The quantity C/'-"^^ is a vector in the 
group space which determines the non-Abelian charge evolution and plays a 
similar role for the charge as J-^ plays for the momentum. The limiting case of 
neutral particles is characterized by J^^ = ~ Q'--^^ = 0. 

2.1 Electrically charged point particles 

In this case the relevant force is the Lorentz force 

r = (2) 

with 

F,u^V,Ak-VkA,, (3) 

where Fik is the Maxwell tensor. The particle under consideration is regarded 
here as a test particle. Moreover, the charge is constant which renders the third 
equation in JQ) irrelevant. 

2.2 Electrically charged spin particles 

According to 20 the evolution of classical relativistic spin particles is governed 
by the equations 

Dp* e ^, I, 

d7 - TO^- ^ 
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and 
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2m 



(5-2) 



(5) 



Here, is the spin four-vector and g is the gyromagnetic ratio. Equation 
((SJ describes the precession of the magnetic moment. It generaUzes earher 
non-relativistic equations by Thomas and Bloch [22, which rely on the 
circumstance that the "expectation value of the vector operator representing 
the "spin" will necessarily follow the same time dependence as one would obtain 
from a classical equation of motion" (cf. HOI)- While the particle momentum 
according to Eq. I^J is independent of the spin vector, the dynamics of the 
latter is coupled (at least via the covariant derivative) to the particle motion. 



2.3 Isospin particles 

Here we have a triplet /^"^^ of scalar fields representing a vector in the three- 
dimensional isospin space, i.e., (A) = (1), (2), (3). This space has an euclidean 
metric G'(^)(s)- The relevant force 



(6) 



where g is the interaction constant, has been obtained by Kerner j^S] and Wong 
pi]. The isospin dynamics is determined by Wong's equation |5J 



where we have used that the structure constants for the SU(2) group coincide 
with the three-dimensional Levi-Civita symbol £''^\b){C)- quantities A^'^'' 
are the vector potentials in terms of which the Yang-Mills field strength tensor 
r},. IS given by 



jk 



The Yang-Mills field equations are 



9£ 



{K){Ly 3 



AT 'A 



(8) 



9- 



(A) 



jk 



(A) 

9£ . (B)(c) 



= 



(9) 



Again we consider the particle motion in a given external field. Wong's equation 
represents a non-Abelian generalization of the equation of motion for electrically 
charged point particles. It can be obtained as the classical limit from quantum 
field theory for the case of sufficiently localized quantum states of the matter 
fields with characteristic length scales much smaller than those associated to 
the Yang-Mills fields [^1^ . 
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2.4 Particles with color charge 

For test particles with color charge the force JT' is given by 

r = /g(^)G(^)(B) , (10) 

m 

with the field strength tensor 

^ife - Vj^^: ^ ^fe^j +.9/ . (K)(L)^j ^fc ' (.11) 

where the / (k)(l) ^'"^ structure constants of the SU(3)group. The field 
equations are 

V.4"'+.9/^.^)(c)A^"^4"T-0. (12) 



9^' 



The quantity Q^"^-* is the color charge with (A) ~ (1) .... (8). Wong's equation 
in this case reads 

(A) 

The structure constants / (b){c) ^"^^ characterized by the commutator relations 

[A(^),A(s)] -2j/(^)(b)(c)A(^\ (14) 

where A(^) are the traceless, hermitian Gell-Mann matrices (see, e.g., |27l I28p . 
In detail we have 

/(1)(2)(3) = 1 , /(4)(5)(8) = /(6)(7)(8) = ^ , 

1 



/(l)(4)(7)-/(2)(4)(6)-/(2)(5)(7) - /(3)(4)(5)— /(3)(6)(7) — /(1)(5)(6)- 2 



(15) 



Below we shall also use the completely symmetric coefficients d(A){B){c} of the 
basic representation which are given by the anti-commutation relations 

A(B)} = ^S(^A){B) + 2d(A)(B)(c)A''^^ , (16) 

where 

C?(l)(4)(6) C'(l)(5)(7) = C^(2)(5)(6) = '^(3)(4)(4) = rf(3)(5)(5) = 
= -'^(2)(4)(7) = -cf(3)(6)(6) = -C?(3)(7)(7) = ^ , 



d{m)(S) - rf(2)(2)(8) - C'(3)(3)(8) - -cJ(8)(8)(8) - -2rf(4)(4)(8) 

1 
3 



2o?(5)(5)(8) = -2rf(6)(6)(8) - -2(i(7)(7)(8) - — . (17) 
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3 Yang-Mills fields with "parallel" potentials 



It is known that for each solution of the general relativistic source free Maxwell 
equations one can construct a set of solutions of the general relativistic Yang- 
Mills equations PHj- Following Gal'tsov ^31 we will refer to the corresponding 
Yang-Mills potentials as "parallel" potentials. The latter are characterized by 

, , g(^)g(s) = l, const. (18) 

Due to the antisymmetry of the structure coefficients the relations (|SJ) and Hll|l 
as well as the equations (O and (|12|l reduce to the linear Maxwell-type forms 
Nevertheless, compared with Maxwell's theory there exists an additional 
degree of freedom, namely the direction of the vector g^"^-* in the group space 
I^U) . Additionally, the structure coefficients f^^(x)(L) different from zero 
which will result in a qualitatively different dynamics. 

3.1 Isospin particles 

In this case the ansatz ()18|l transforms the first of equations ^ with (0 into 

which has the structure of the equations of motion of a particle with charge 
e = gl^^'^qiA) under the infiuence of the Lorentz force Analogously, one can 
rewrite equation Q for the isospin evolution, 

= -n e^'Uc^^h^-^ , n . £a.,^ . (20) 

Because of the antisymmetry of the Levi-Civita symbols the equations H2()|l ad- 
mit a quadratic integral of motion /'^'^^/(^^ = const, which is a Casimir invariant 
|22| . normalizable to /^"^'/(^) = 1. In addition, we obtain from l)2()|l 

g(^)^— = 0, /(^)g(^) = / = const . (21) 

Using the standard definition 

{A) 



I . ^ 



of the vector product, Eq. H2UI) may be written as an equation for the precession 
off. 



I , n 



(23) 



The "longitudinal" component l''^^q{A) i the projection of the dynamical variable 
l'^^^ on the "rotation axis" q'^^\ remains constant according to H21|) . 
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3.2 Colored particles 



Similar to the previous isospin case the color charge evolution equation 
admits the existence of a quadratic integral of motion Q^^'^Q^a) — const which 
is the first Casimir invariant [221 ■ The condition (|18|l of parallelism in the color 
space provides a second integral of motion — const, as well. However, 

the corresponding 8-dimensional evolution equations 



H 



(A) ^ j,{A) 



(24) 



are more complicated than the equation H23|l for the isospin precession. Different 
from the SU{2) model there exists a second Casimir invariant 



(25) 



where d 



(A){B){C) s-re the totally symmetric group coefficients IjlTf) of the given 
ion 

' (g(«))^ + X/3Q(«) (Q(1))2 + (Q(2))2 + (Q(3))2 



basic representation of the SU (3) group. In detail it reads 
1 



+3g(2) 



(g(^))^+(g(^))^+(g(«))^+(g(^))2 

3 



~3g 



(1) 



)(*)g(«)+g(^)g(^) 



-g(^)g(^)+g(^'g(«) 



(3) 



(Q(4))2^(Q(5))2_(g(6))2_(Q(7))2 



(26) 



4 Gravitational wave background 

Our aim in this paper is to study the general dynamics outlined so far in the field 
of a plane- fronted GW with parallel rays (PP-wave). We assume the latter to 
be an exact solution of Einstein's vacuum field equations with a five-parametric 
group of isometrics G5, including a covariantly constant null Killing vector (KV) 
|32| . Gravitational waves are usually described either in Fermi coordinates or in 
the transverse-traceless (TT)-gauge. For the merits of each of these choices and 
for issues of gauge-invariance in the linearized theory see, e.g., [S^. In order 
to establish a comprehensive picture we start by sketching our basic setting for 
both cases. For computational ease most of the analysis will then be done in 
TT-coordinates. 

4.1 PP-wave in Fermi coordinates 

The corresponding line element 

ds^ = 2dudv - dy^ ~ dz^ - 2H{u, y, z)du'^ , u = , v = , (27) 
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contains a harmonic function 7i, obeying 

which is quadratic in y and z for a symmetry group. Exphcitly, we have 

2H(u, y, z) = A{u) {y^ - z^) + 2B{u) yz , (29) 

where A{u) and B{u) are arbitrary functions of the retarded time u. We may 
define a periodic GW by assuming the variables A and B to be periodic functions 
of u. All non-vanishing components of the Riemann tensor, Rzuzu = —Ryuyu = 
A{u) and Ryuzu = —B{u), are periodic for this case. The GW metric in Fermi 
coordinates is non-singular for arbitrary retarded times because det{gik) = —1 ^ 
0. However, the weak field approximation max|(/ifc| << 1 is correct only close 
to y = 0, z = 0. 



4.2 PP-wave in TT-gauge 

The line element in TT-gauge has the form 

ds^ = 2duAv - [cosh27 {e^f^{dx^f + e-^^{dx^f) + 2 six\h2-i<lx^ dx^] , 

(30) 

where u = (t — x^) 1^/2 and v = {t + a;^) /-\/2 are the retarded and the advanced 
times, respectively. For this metric the three KVs which form an Abelian sub- 
group of G5 are 

C^.)=5;, e^2)=^2, ^(Z)=^\- (31) 

The KV ^^^j is a covariantly constant null vector, orthogonal to ^^2) ^^"^ C(3)- 
The functions /3(m) and ^(u) are arbitrary. We shall focus here on the case 
of periodic functions /3(w) and 7(u). This definition of periodicity docs not, in 
general, coincide with the definition in Fermi coordinates given above. However, 
both concepts of a periodic GW have the same weak- field limit, which in the 
TT-gauge is characterized by L = 1, |/3(m)| << 1 and |7(u)| << 1. Generally, 
the function L(u) satisfies the Einstein equation 

L + L ((/3)2 cosh^ 27 + (7)2) = , (32) 

where a dot denotes the derivative with respect to the retarded time u. We 
assume the hypersurface m = to be the leading front of the GW with 

/3(0) = 7(0) = , i(0) = l, /3(0) - 7(0) = i(0) = . (33) 

For the special case 7(7/) = (equivalent to B (u) = 0), corresponding to only 
one polarization direction, the transformation relations between Fermi and TT 
coordinates are 

u ^ u, v^v^-^{{x'f{L'e^f')■^-{x'■f{I?e-^^)■\ , 

y = Le^•x^ z = Le-^-x^, A(m) = /3 + 2/3y . (34) 
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The last formula clarifies the relation between the different periodicity defini- 
tions given above. For the physically relevant situation where the background 
factor L changes only slowly compared with the change of the wave factor (3 [cf. 
Ref. the last term in the formula for A{u) in Eq. (|34|l can be neglected. 

Then we have A{u) = (3 and both periodicity definitions coincide. 



5 Particle dynamics: Solutions 

We are interested here in the particle dynamics in given external gravitational 
and Yang-Mills fields. The restriction H18|l to Yang-Mills fields with "parallel 
potentials" simplifies the Kerner-Wong equations since the quantities 
and Q'^^^q(A) remain constant. We introduce the cumulative symbol a for either 
e, or gl^'^'^q^A)! or gQ^^\(A)i which allows us to write the equation of motion 
with cither or ®, or (|10|) in the unified form 

^--F\P\ ^ = ^. (35) 
Dt m (It m 

The orthogonality of the force to the particle momentum corresponds to the 
quadratic integral 

9'''p^Pk = m2 . (36) 

One may solve this relation for one of the components of the momentum. In 
Fermi coordinates we have 

Pu^:^ [m^ +pI+pI- 2H{u, y, z) pi] (37) 
The analogous formula in TT coordinates is 

Pu^T^[m^-g''^{u)pc.Pf3] , (38) 

ZPy 

where greek indices run from 2 to 3. Both in Fermi and in TT coordinates the 
covariantly constant null KV has the form ^^-^ = 61 and ^^^^ = 61, respectively. 
Below we shall restrict ourselves to fields which satisfy 

q.fF^k ^ q,-,F,k = . (39) 

This implies that the quantities C(^^^Pi ^-nd ^(^•jPi are integrals of motion (see 
e.g. ESI), 

C(v]P^ = C(v)Pi ^ Pv ^ Cv = const . (40) 
Using the general relationship 

m—=p^=p''^p,^C^, (41) 
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one can reparametrize the remaining equations for ^ 0, by means of the 
hnear formula (notice that u — u [cf. Eq. (|34|l ] holds also in the general case) 

771 

T = Tn + —u . (42) 

We start our solution procedure by first recalling the reference case of neutral 
particles. 



5.1 Neutral particles 

For a vanishing generalized charge a the equations of motion in TT coordinates 
are immediately integrated. The result is 

pM = a , P2{u) = C2 , pz{u) ^C3,pu = -^ [m^-g^^iu) C^^Cp] , (43) 

where C2 and C3 are constants. A particle moving in direction of the GW 
propagation before the infall of the latter, i.e., P2 = Pa = will not change its 
direction. If we additionally have C„ = the particle is at rest both before 
and after the GW infall, since pi{u) = P2{u) — p^{u) ~ and po = mc. An 
observer at rest, characterized by a four- velocity V"^ — -^(<5,\ + S\,) = S'q, would 

measure the (invariant) particle energy £ = p^Vk = '^{Pu +Pv)- For neutral 
particles Pv and Pu are given in Eq. H43|) . Corresponding expressions for charged 
particles will be obtained below. 

In Fermi coordinates the situation is as follows. The system (|35() reduces to 
the set of equations 

y = -C^^py , i = -c~ V , 
C-^Py = -A{u)y-B{u)z, C'-'^p, ^ A{u) z - B{u) y , (44) 

for y, z, py, and p^. If the latter set of quantities is known, the component Pu 
follows via H37|l . The quantity v{u) may be found by solving the equation 

* = C-Vu + 27i(u,y,z) . (45) 

The set (|44|) is a linear, homogeneous first-order system of differential equations. 
For a GW field with polarization B{u) = and with A depending on the 
retarded time via the dimensionless variable ku, it can be written as 

y" - Aiku) y = , z" + A{ku) z = . (46) 

Differentiation with respect to ku is denoted by a prime and A (ku) = A (u) /fc^. 
For A{u) = 0, i.e. in the absence of a gravitational field, y{u) and z{u) are linear 
functions of the retarded time (and of the affine parameter r) , and the particle 
has constant momentum. For a periodic function A(u) the equations H4t)|) are 
of the type of Hill's equation PIEIIHI- The solutions are Hill functions. For 
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a dependence A{ku) ~ 5 + ecos{ku) where S and e are constants, Eqs. 
reduce to Mathieu equations which have solutions of the type 131 |S| 



y (X exp [fj,ku\4> (ku) + cxp [—iiku\ijj (ku) , 



(47) 



where (j) and ^ are periodic functions with the period of A{ku), i.e. in the 
present case, (j){ku + 2TT) — (l>{ku) and ip {ku + 2tt) = il}{ku). The solutions 
consist of products of an exponential function and a periodic function of period 
277. The characteristic exponent /i is a complex constant. For Re {^) = the 
solution is stable. In general, it is not periodic again but it is oscillating (T', 
p. 115). For Re {(i) ^ either the first or the second exponential function in H47|l 
is unbounded and the solution is unstable. The stability properties of Mathieu's 
equations are well known in the literature and may be visualized by stability 
regions in an 5 - e diagram (see Eq.(4.1) and Fig. 5.1). Let's consider the 
stability region which is closest to the origin 5 = e = 0. For small positive values 
of 5 and e the boundary of this region is determined by the line 5 = ^ — ^e. 
Apphed to the case A{u) = AQCOs{ku), i.e. (5 = and e — A^/k'^, this means 
stable solutions for e = ^ < 1/2. Since Aq — (3 (0) — (i^k"^, we have also 
/3o < 1/2. Under this condition neutral particles are parametrically oscillating 
in Fermi coordinates. 

While equations of the type of Mathieu's equation and questions of stabil- 
ity will play an essential role in the following investigations of the dynamics of 
charged particles (see subsection 15.2.21 below) . it is obvious that the description 
for neutral particles is more involved in Fermi coordinates. Therefore, for com- 
putational ease and in order to separate charge effects from the neutral particle 
motion we shall perform the following analysis in TT coordinates. 

5.2 Charged particles 

To obtain exactly solvable models for the particle motion we resort to simple 

(A) 

field configurations F^f, . For particles with electric charge we focus on the 
motion of the latter in a constant homogeneous magnetic field Hq orthogonal 
to the GW front plane, which corresponds to a Maxwell tensor 



A corresponding generalization for non-Abelian fields with parallel potentials 
according to ((TH)l is 



Both and (03 satisfy with . This constitutes a model in which the 
gravitational wave and the fields H48|l or 149|l are given, external fields which are 
independent of each other. It may provide the basis of a perturbative treatment 
with respect to the GW amplitude within a linearized theory. It is worth men- 
tioning, that the expressions (|48|l and H49() are also solutions of the Maxwell- 
and Yang-Mills equations, respectively, on the background of the exact GW 




(48) 



f'^^ = <j(^)Af {5^^6l - 5pl) , M = const . 



(49) 
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(or (|77|l ). This allows a study of the corresponding field dynamics on a GW 
background, which, however, is not the purpose of the present paper. 
For the equations of motion in TT coordinates we obtain 

dp2 Ma .32 , 33 ^ dps Ma .32 , 23 \ ,rn\ 

-5— = [9 P2+9 P3) , -5— = (.g P2+9 P3) ■ 50) 

dr TO dr TO 

Equivalent second-order equations are 

d^P2 , n /,.xdP2 , TTW..X„ n 1/^1,1 23 



^ + i?3(^^)^ + Vl^3MP3 = , p2 = -L(^_^p3_g23p3) , (52) 

where we have introduced the notations 

i?2(u) = 















g'^u) 






522(u) 







and 



Rsiu) = _il^^2 -+/?-7tanh(27) , (53) 



^2(«) = ^ + ^f/3sinh(27) ^ 



with 



L'^ \ ' " cosh(27) 

^3(.) = + f (/^-M27) + ^) , (54) 

n = — ^ = const . (55) 

The substitution /3 — > — /3 converts i^s into i?2 and vice versa, while W3 is 
obtained from W2 by /? ^ — /3 and simultaneously 11 — > —11. 
By the substitution 

p„ = Z,(M)exp|-i^ i?a(CK| (56) 

the equations H51|l and H52|l may be transformed into the Hill equations 

Z^ + Fc,iu)Z^^O , (57) 

where 

F^^W^-^-^. (58) 
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The detailed form of the relations (|56|) and H58|l is 



P2=^2(u)v/cosh(27) — , p3=^3HVcosh(2^ — , p„(0)=Z„(0) = 



(59) 



and 



2n 



+ /3 + 7 tanh(27) 



^^^^^^ + (/3)2(sinh(27))2 + 2^/3 + 2^7tanh(27)-27/3tanh(2^ , (60) 



Hi? 



n2 2U 



(i sinh(27) 



7 



cosh(27) 



/3 + 7tanh(27) 



■ ( x.(o ^^2+(/3)'(sinh(27))2-2^^/3 + 2^^7tanh(27) + 27/3tanh(27) , (61) 
(cosh(27))"' L L 

respectively. In a next step we have to solve Hill's equations (|57|l . 
5.2.1 General solution of Hill's equation 

The structure of the solutions of the linear, second-order differential equations 

(EH) is ElEllTllHlEl 



The functions Ha (u) satisfy the initial conditions 

i/2(0)-l, ff2(0)-0, H3{0) = 0, H3{0)^1 



(62) 



(63) 



and represent the fundamental solutions of Hill's equation (I57|l with unitary 
Wronsky determinant. For Z:}{u) we have 



Z^iu) = C3 H* + C2H , 



(64) 



where 

m 



IP cosh 27 

L2 



H2 + H2 ( 7tanh27 + - J - sinh27 



L 



H 



H3 + H3 I 7tanh27 + - ^ J2 sinh27 



cosh 27 

(0) = i?2 (0) - , h; (0) - i/3 (0) = 1 



(65) 



(66) 



(67) 



In the absence of gravitational radiation, i.e. for /3 = 7 = 0, L = 1 , the functions 
Fa in equation lfS5|l reduce to 



F2{u) = Fsiu) = const = 



(68) 
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Equation 157|l then describes harmonic oscillations with 

H2 = i?3* = COS Uu, H3 = iJ* = ^ sin Uu . (69) 

Since with (|42|l and H55|l we have IIu —> Qht where fin = is the Lar- 
mor frequency, we recover the corresponding particle rotation in flat spacetime. 
Generally, the functions H2, H^, H2 and cannot be written in terms of 
elementary functions but are given as series representations. In the following 
subsection we shall be interested in expressions for F2 and F3 for which the Hill 
equations H57() specify to Mathieu equations. Then the functions Ha can be 
expanded in powers of the GW amplitude, where the zeroth order is given by 

(inni). 

The analysis so far may be summarized by writing the solution of the equa- 
tions of motion (|51|l and H52I) in the compact and elegant matrix form 

where 

^ ^/J2^ f \ f H2{u) -UHs{u) 

L [ e-^ )■[ nHiiu) H^{u) 

The set of equations H7l)|) and (|71|l represents the general solution for the mo- 
mentum of the charged particle in the GW field H3U|I and the Yang-Mills field 
(133). 




5.2.2 A simple model 

Now we apply the general formalism to a "sandwich" GW (see, e.g., [33]) with 
polarization 7 = 0. Let (3 be different from zero during a finite retarded time 
interval T, i.e., /3 = for u < and u > T. Within the interval < u < T we 
assume /3 to be periodic according to 

/3(u) =/3o(l-cos(fcu)) , (0<u<T). (72) 

This implies /3(0) = /3(27r/fc) = and ${0) = /3(27r/fc) = 0. Furthermore, 
the time scale T is assumed to be small compared with the scale on which the 
background factor L changes |33]. Under these conditions we may neglect the 
L/L terms in (|6U|I and (|51|l and use the latter expressions with L = 1. For this 
situation the potentials F2 and F3 reduce to 

F2 = Pok^ cos (ku) = - i?2^2„ , 

F3 = tf-/3ofc2cos(M =n2--Rl3« (73) 

in the interval < u < T. Replacing now the variable u by ku and denoting 
the derivative with respect to ku again by a prime, the equations (|57|l with H73|l 
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specify to 



Z'i + + /3o cos (fcu)^ Z2 = , 



^3 + - /9o cos (fcu)^ Z3 = . (74) 

Both Z2 and Z3 obey Mathieu equations. In the absence of the GW, i.e. for 
/3 = 0, we have Fa{ku) — Il^/fc^ = const (here we have used the redefinition 
Fa (ku) = Fa (u) /k'^ ) and the equations of motion reduce to harmonic oscillator 
equations with solutions 

Z2=P2= C2 cos (nu) - C3 sin (nu) (75) 

and 

Z3 = p3 = C3 COS (Jlu) + C2 sin (Hu) . (76) 

Replacing here 11 and u according to (|55|l and H42I) . we find a particle rotation 
in the x^Ox^ plane with the angular velocity flM = which is, of course, 
the analogue of the Larmor frequency. Immediately after the wave front, i.e. at 
u = 0+0, we have 

F2{0) = ^ + Po, F3{0) = ^-f3o. (77) 

The jump Pok'^ of the curvature tensor at the front makes the evolutions of the p2 
and p3 components different. They begin to oscillate with different frequencies 
and the particle trajectory is no longer circular. Corresponding features hold 
for the second polarization f3 = and 7 7^ 0, for which F2 and F3 differ in the 
term hnear in 11 [cf. Eqs. l|Bn)l and 

The general solutions of Eq. ifZ^ are of the type of "cosine elliptic" and 
"sine elliptic" functions (see [3 El)- As already mentioned, the latter may be 
expanded in powers of /3o with the zeroth-order terms 1)75(1 and ((76|l . 

Eqs. (|74|l are of the same type as Eqs. (|46|l . With the identifications S 
iP/k^ and e f3o, the stability discussion mentioned in subsection 15 . II mav be 
apphed here as well (see 0, Eq.(4.1) and Fig. 5.1). Depending on the parameter 
combinations the solutions may be stable or unstable. Within the stable regions 
the functions Z2 and Z^ are parametrically oscillating which, according to H59|l 
implies a corresponding behavior of the particle momenta. The regions of stable 
solutions are characterized by stability zones in a (Il/fc)^ x (3q plane which are 
connected together at the points (Il/k)'^ ~ ri^ f3o ~ 0, where n is an integer 

El 

The analysis of the neutral particle motion in subsection 15 . II corresponds to 
the case 11 71 = 0. For 11 7^ the relevant values for the transition points 

are n ~ 1,2, For n = 1 we have H/k — 1/2. These points on the axis 

Po = (which corresponds to the absence of the GW) are the only transition 
points between stable regions which also belong to the stable region. All other 
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boundary points of the stable regions are unstable points. Consequently, any 
deviation from /3o = 0, i.e. even a GW with arbitrary weak amplitude f3o, 
induces an instability in these critical points. In particular, this is true for 
the point H/k = 1/2 (see [2], Fig. 5.1). This demonstrates that parametric 
instabilities are a generic phenomenon for the motion of particles in all the 
cases considered here. While we have obtained this result in TT coordinates, 
the transformations (|34|l allow us to find the corresponding particle momenta 
in Fermi coordinates as well. The relevant transformations are 



Pv=Pv=Cv, py=p2{u)- 



L 



P\y{u); 



P3{u)- 



L 



./3 z(w) 



(78) 



where 



/(w) = Le^ 



z{u) = Le 



y{0) 



dCP2{0L-\0e-^P(^^ 



m- / dCpsiOL-^Oe^^^'^^ 



(79) 



It is interesting to realize that there are astrophysical situation for which the 
existence of such kind of instabilities might be relevant. This can be seen with 
the help of the following order-of-magnitude estimates. Eq. (|55|l may be written 
as n = Q.M'm/Ci, with = Majm. For the electromagnetic case one has 
(7 = 6 and M = Hq. The interstellar magnetic field is of the order 3 -;- 6 • 10~^Oe 
IS]. The integral of motion is equal to C„ = (v^m^ +p2(o) - pi(0))/v^. 
For nonrelativistic particles oc m/\/2 and 11 cx ojrV^- (The coefficient 
disappears if we use the natural parameter r instead of retarded time u). Using 
the estimate cx lO'^ • i?o(Hz) we find ujh « 10^ 10^ Hz. This is well 

within the typical range 1 — lO'^Hz for the frequency /c of a GW, generated by 
rapidly rotating neutron stars (pulsars). Thus, for non-relativistic particles 11 
may be of the same order as the GW frequency k. The situation is different for 
ultrarelativistic particles. Since for particles that move in propagation direction 
of the GW (p^(0) > 0), we have ^ for m ^ 0. Consequently, the quantity 
n becomes very large [cf. Eq. (|55|l ]. such that 11 >> /c. However, if the particles 
move in the opposite direction, i.e. p^(0) < 0, any value of Cy is possible. In 
particular, H is not excluded to be in the range of about lO^'^Hz which is typical 
for infra-low-frequency GW from relativistic compact binaries. 

Since an ensemble of particles will generally not be characterized by a single 
value of Cy but by a distribution, the quantity may play the role of a tuning 
parameter in the following sense. Let us associate a mean value (C^) for the 
system as a whole and let the system be outside the resonance uiH/k — n/2 for 
this (Cy), but not very far from it. Since H = ujHm/Cy, we will very likely find 
a particular particle with a specific Cy such that for this particle H/fc — n/2 
exactly. Consequently, certain particles of the ensemble might be resonant under 
these conditions. 
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6 Evolution of the spin four-vector 



In this section we focus on the solution of Eq. (O in external gravitational 
and magnetic fields. The equations Q and (0 admit two integrals of motion, 
namely piS^ = and SiS^ = — = const (20j. In the present context this 
amounts to 

PvSu + PuS, +Pc,S'' = 0, 25„5„ + Sc.S" = , (80) 
which may be used to eliminate the components Sy and Su according to 
1 



2pv 



(81) 



and 



Sv — r, 

2Pu 



PaS" ± J{po.S'-f + (to2 _ p^p'-){El + Sc.S'-) 



(82) 



respectively. After the reparametrization 142(1 the remaining equations are 



Si: 



= i^g^'gpa ( s,-p,^ ) + — 



gF^k + ^^-^PcFmS'^p' 



(83) 



Furthermore, the equation for Sy is 



dSy _ e{g - 2) ^ ^ , 



du 



2m2 



-FuS^p' 



(84) 



In the following we shall restrict ourselves to the exactly integrable case g = 2. 
Under this condition we find from 184|l that 



Sy — const = Ey 



(85) 



After the substitution 



Sa — -^a ~t" ^ Pa 

L^y 



(86) 



where Pa is assumed to be given by Eq. (|70() in terms of Hill (or Mathieu) 
functions, we obtain a homogeneous equation for the new variable X^, 



— ] - 77.9 9paX„ + —taag 

au 2 Ly 
It is convenient to write this equation in the matrix form 



(87) 



du 



A+-^B|.X 



(88) 



where X is a column vector with elements X2 and X^. The two-dimensional 
matrices A and B have the structures 



^ +cosh2(27)/3 
3-2/3(^ + sinh(27) cosh(27)/3) 



62/3(7 - sinh(27) cosh(27)/3) 



cosh'^(27)^ 



, (89) 
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and 

/ sinh27 -e^/^ cosh 27 \ , , 

^- 1^ 6-2/3 cosh 27 _sinh27 J' 

respectively. Equations (|50|l may be written in matrix form as wcU: 

^O'.ilB.f . (91) 

For the derivative of the combination + P3 obtain 



^ {pI +pI)=^ M - pI) sinh 27 - 2p2P3 cosh 27 sinh 2/3] . (92) 



In general, the right-hand side of this equation is different from zero, i.e., the 
particle motion is no longer circular in the GW field. 

Eq. (|88|l may be further simplified by changing to a new variable Y, defined 

by 

X = T • Y , (93) 

where T is supposed to satisfy the differential equation 

-^T = A • T . (94) 

du 

This procedure [cf. |35[ I36| ] removes the A term in Eq. (|88|l and gives rise to 
the equation 

for Y with 

B = T"^ • B • T . (96) 

By direct calculation one checks that the matrix 

T = L-( '^"t ^ ^ V f '"''f " '"/"^ ) > (97) 

\ 6 / V smh7 cosh 7 J y smip costp J 

where 

u 

V' = y /3sinh27d'u , (98) 


satisfies the equation H94(l . The determinants of each of the three two-dimensional 
matrices in H97|l are equal to one. In the absence of the GW field all of them 
are identical to I, i.e., 

T(0) = I^ ( J J ) . (99) 

The structure of the third matrix on the right-hand side of Eq. (|97|l suggests 
the interpretation as a gravitationally induced rotation with phase ipi''^) s-nd 
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frequency ip{u). For either of the polarizations 7 = or /3 = 0, however, we 
have = and the third matrix reduces to I. 

Direct calculation of the matrix B in (|96|l with the help of the expressions 
H9U|) . H97|) and (|98|l yields the surprisingly simple result 

B = • B • T = J ) ■ ^^^^^ 

It is remarkable, that the matrix B, different from B, does not depend on 
retarded time. This property allows us to find the solution of equation 195|) in 
terms of elementary functions as 

¥3(1*) J V ^3(0) / V sm$(M) cos$(u) y 



where 



The combination 



u 



Y,'{u) + Y,\u) = ^2^(0) + y/(0) (103) 

is preserved, i.e. the dynamics of Y represents a rotation in the x^Ox^ plane. 
The functions iS'2(u) and 5*3 (u) in H86|l ) can now be expressed in terms of the 
three matrices H(u), T{u) and R(m), given by the expressions (|7T|) . IpTfl and 
HlOljl . respectively: 

IS ) =-^<«> ■ ■ ( t(o) ) +1; I «<")-T<"> ■ J'l")! ( § 

(104) 

While the matrices T(u) and R(m) are constructed out of elementary functions, 
the matrix H(it), according to Eq. (|71(l . consists of Hill functions, which for 
the special case of subsection 15.2.21 reduce to Mathieu functions. The latter, in 
turn, can be expressed via "cosine elliptic" and "sine elliptic" functions (see the 
discussion following Eqs. 1)77(1 ^. All these functions are assumed to be known 
here. In the absence of the GW, 

T(.) ^ I, H(.) = R(.) ^ Ro(r) = ( ^^^^^ "jl^i^^J ) , (105) 



sinr2/fT cosHht 
and we recover the standard flat spacetime rotation of the spin particle. 

Equation 1] 104(1 represents the general solution for the spin dynamics in the GW 
field (|30() and the magnetic field H48|l . The structure of the solution 1(1 04() allows 
us to interpret the spin dynamics as composed of three separate contributions. 
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characterized by the matrices R(w), T(m) and H(u). The matrix R represents 
a Larmor type precession with the frequency j^^- As already mentioned, 
the matrix T describes a gravitationally induced rotation with phase ij^iu) and 
frequency ip{u), and finally, the matrix H accounts for the coupling of the 
particle motion (|7()|l to the spin dynamics. 



7 Isospin evolution 

In this section as well as in the subsequent one we discuss the dynamics of 
non-Abelian charges under the influence of external gravitational and Yang- 
Mills fields. Let us consider Eq. H2U|) for the SU{2) symmetry group. Since the 
structure constants for this group coincide with the Levi-Civita symbol, all three 
directions in the isospin space are equivalent. With the choice l'^"^^ = /(a)'?^'^'* 
we obtain the following equations for the isospin evolution: 

^ = n.ii^K ^ = -^-/«. (107) 

dr dr 

The precession frequency fl in Eq. pOfl is calculated on the particle worldline 
with from H7()|l and the potential 



(108) 



corresponding to the constant solution — M . Differentiating the expression 
pU8|) . we recover the field strength H49I) . The frequency il(u) in H2U|) is given by 

n{u) ^ ^{A2P^ + Asp^) . (109) 
m 

Here, the arbitrary constant was chosen such that 17(u = 0) = 0. In order to 
find the terms a;^(u)— x^(0) and x^{u)—x^{0) which are needed in (|108|1 . we have 
to integrate the second equation in H35|l. The formal solution is 

x"(u)-x"(0) = d^g'^^iOppiO , (110) 

which provides us with 

"^"^ ^ 2^ {/"^^ [p'(OpHu)'PHOp'i^)] - ^ [C2p'iu) - Csp^u)]^ . 

(Ill) 

Again we assume here p^(u) and p^{u) to be known, i.e., the particle dynamics 
is considered to be solved [cf, Eqs. (|70|) and (|71|l ]. The solution of the system 
<|1U7|) then becomes 



= /cos\E'(m) , = -/sin*(w) , *(w) = ^-(0) + ^ / 



D,{u)du . 
(112) 
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The function ^{u) plays the role of the (generally M-dependent) phase of the 
isospin precession in the external Yang-Mills field |2Z|- The set of equations 
H112|) with provides a complete description for the isospin dynamics under 

the influence of the GW (jSOJl and the Yang-Mills field 



8 Color charge evolution 

The SU(3) case may be studied along similar lines although it is technically 
more extended since more degrees of freedom are involved. As a result, we shall 
find a richer dynamical structure than in the SU(2) case. With the ansatz H18|l 
and the expressions (|15|l . Eqs. H24|l for the color charge dynamics become 



dQW 

dr 
dQ(^) 

dr 

dg(^) 

dr 
dQW 
dr 



dr 



dQ(6) 
dr 



dQ(^) 
dr 



dQW 
dr 



-n 

n 
n 



(g(^)g(3)-5(3)Q(2))+l(^(4)g(7)_^(7)Q(4))_l(g(5)Q(6)_^(6)Q(5)) 
(g(3)g(l)_g(l)Q(3))+l(^(4)Q(6)_^(6)Q(4))^l(g(5)Q(7)_^(7)Q(5)) 

(g(i)g(^)-9(2)g«)+i(gWg(5)-g(5)gW)-i(<z(6)g(7)_^(7)Q(6)) 



(q(r)Q(^)-qWQi^)) + {q(^)Q 



(2)_„(2)r)(6) 



n 

2 



(g(l)g(«'-g(6)g(l))^(^(7)Q(2)_^(2)Q(7))^(^(3)Q(4)_^(4)Q(3)) 



n 

2 



((7(5)Q(l)_5(l)Q(5)) + (g(2)Q(4)_^(4)g(2))^(^(3)g(7)_^(7)g(3); 

[g(^)g(«)-g(«)g(7)] , 



n 

2 

^/3^2 



(<7(l)g(4)-g(4)g(l)) + (9(2)Q(5)_g(5)Q(2)) + (5(6)Q(3)_g(3)Q(6)) 

[gWg(6)-9(6)Q(8)] , 

(^(4) g(5) Q(4) ) + (g(6) Q(7) g(6) ) 



(113) 



The space of color charges may be split into three different subspaces, which 
correspond to the structures of the SU{2), SU{2) x U{1) and U{1) subgroups of 
the total group SU{3) (see, e.g. US]). In the following subsections we consider 
the vector g'"*' to lie in the first, second, and third subspaces, respectively. 
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8.1 First special case 

Let the vector q^^^^ have only the three non-zero components q'^^\ g*^^^ and q^^\ 
It is then evident that 

(g(i))2 + (g(^))2 + (g(3)f = const , (iw) 

(Q(^))' + [Q^'^r + {Q'^'^f + {Q'^'^f = const , (115) 

and 

g(*) = const . (116) 

For the color charges Q'^^\ Q^^\ and Q''^\ which correspond to a SU{2) sub- 
group of the total SU (3) group, the combination pi4|l is preserved. A similar 
relation holds for the set Q('^\ Q''^\ Q^^\ and Q^'^\ while Q^^'> is separately 
conserved. Relations ()114ll - H116ll are also obtained for the case that the only 
non- vanishing components are q^^\ q^^\ <l^^\ and q^'^\ as well as for the choice 
— g(2) = = q(7) _ Q g^^^j ^(8) _^ YjQi us now further specify to the case 

q^^) — 5'"f^fy Then the system (|113ll takes the form 

dr dr 



= 'no.'. , 'jh. = _ 'no... , 

At 2 At 2 

The color charge g*^^-* remains constant because it is the projection Q'^^^q(A) 
of the vector g^"*^ on the given preferred direction q(A)- The charge g^*-* does 
not evolve, because for such a q^A) the antisymmetric tensor H(a)(B) in H24|l 
does not contain a non- vanishing component with (A) = (8). The equations for 
g^^-*, ... g^''-' split into f/iree two-dimensional subsystems with the pairs g*-^^ 
and g(3), g(4) and gC^), g(5) and g(^). The evolution of the first pair (g^^) 
and g*-^-*) corresponds to a precession in the group space with the frequency f2. 
It has a solution of the type (|112|l . The dynamics of the pairs Q^'^\ g^^-* and 
g(^), g(^) is a precession with the frequency Q/2. 

8.2 Second special case: g*^"^^ = 6^ff^ 

Now we assume g'^'^^ to lie in the second subspace. As an example we consider the 
case q^^'^ — S'^^fy For this choice the set of equations (|113|l can be transformed 
into 

dgw 

dr 
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= i§i^ = ><^>. (118) 

dr 2 dr 2 

The quantities Q^^^ and — + ^Q^ remain constant, the projections Q^^^ 
and Q* = ^{Q^^^ + v^Q^^^) precess with the frequency Q, the pairs Q^'^\ Q^^) 
and Q^'^\ Q^^^ precess with Q/2. Furthermore, one has 

{Q^'^f + {Q^'^f = const , (Q(2))2 + {Q^'^f = const , 

(Q^^^f + (Q^'^f + {Q^^^f = const . (119) 

8.3 Third special case: g(^) = d'^^^ 

Finally, let q^"^^ be parallel to the basis vector in the U{1) subspace, i.e., q^^^ = 
. Here we obtain 

dg« ,^ dO(2) dQ^_ dQ(8)_ 



dr dr dr dr 



as well as 



The pairs Q^^\ Q(^) and Q^^\ Q^^) precess with the frequency ^O, while Q^^), 
(5(2), Q(3), and are constant. 



8.4 Remarks on the general case 

In the general case one expects the color vector Q^^^ to rotate in the hypersur- 
face orthogonal to q^^^ in the group space. This is illustrated by the following 
analogy. Let us consider the standard decomposition of the Maxwell tensor with 
respect to a four-velocity vector of an arbitrary observer, 

Fik = EiVk - EuV, - e,kjiH^V^ , (121) 

where E"^ and H'^ are the corresponding four-vectors for the electric and mag- 
netic fields, respectively. For FikV^ = the comoving observer experiences 
a magnetic field only which, via the Lorentz force generates a spatial particle 
rotation, i.e. a rotation in the hypersurface orthogonal to V^. 
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In the present case the vector q^"^^ plays the role of V^. Instead of Maxwell's 
tensor we have to consider [cf. Eq. 1241) ] the antisymmetric tensor = 

f{A){B){C)'l^^^ ■ Since -ff(A)(c)'i''''^'' — Oj this represents a group space analogue 
to the previous case of a pure magnetic field in spacetime. Consequently, the 
color vector rotates in the group space analogously to the momentum vector of 
a charged particle in a pure magnetic field. 

9 Conclusions 

Charged particles in electromagnetic fields are known to be parametrically in- 
fluenced by gravitational waves. Typical phenomena are parametric resonances 
and parametric oscillations. In the present paper we have generalized and ex- 
tended work in this field to include classical spin particles and particles with 
non-Abelian charges in specific Yang-Mills fields. Moreover, no weak field 
approximation was used for the GW. The electrically charged spin particle 
was described by the Bargmann-Michel-Telegdi equations. For the dynamics 
of the non-Abelian charges we used Wong's equations for the isospin (SU(2)- 
symmetry) and for the color charges (SU(3)-symmetry). We derived exact gen- 
eral solutions for the parametric influence of the GW on the particle motion 
in each of the mentioned cases, including the dynamics of spin, isospin, and 
color charge. For the case of a special sandwich GW the particle dynamics was 
reduced to a set of Mathieu equations. Using well-known stability properties of 
the latter we found that parametric instabilities are a generic phenomenon for 
such kind of particle motion. Since spin, isospin and color charge are coupled 
to this motion, their dynamics is affected correspondingly. The spin dynamics 
was shown to be composed of three elements, namely a gravitationally modified 
Larmor precession, a part due to the coupling to the particle motion, and a 
pure gravitational part. The vectors of isospin and color charge carry out grav- 
itationally influenced precession motions in their group spaces which we have 
classified for several cases. 
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